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Abstract 

Heavy quark production in 7 TeV centre-of-mass energy pp collisions at the LHC 
is not necessarily flavour symmetric. The production asymmetry, A-p, between 
and mesons is studied using the (/>7r^ decay mode in a data sample of 
1.0 fb~^ collected with the LHCb detector. The difference between vr"*" and vr" 
detection efficiencies is determined using the ratios of fully reconstructed to partially 
reconstructed D*^ decays. The overall production asymmetry in the rapidity 
region 2.0 to 4.5 with transverse momentum larger than 2 GeV is measured to be 

= (—0.33 lb 0.22 lb 0.10)%. This result can constrain models of heavy ffavour 
production. 
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1 Introduction 



Production of charm and bottom hadrons at the LHC in 7 TeV pp colhsions is quite prohfic. 
The bottom cross-section in the pseudorapidity region between 2 and 6 is about 80 M-b [I] , 
and the charm cross-section is about 30 times higher p]. In pp colhsions the production 
rates of charm and anti-charm particles need not be the same. While production diagrams 
are flavour symmetric, the hadronization process may prefer antiparticles to particles or 
vice versa. Figure [T] gives an example of cc production via gluon fusion. If the quarks 
that contribute to charm meson production are created in an independent fragmentation 
process, equal numbers of D and D will be produced. On the other hand, if they combine 
with valence quarks in beam protons, the c-quark can form a meson, while the c-quark 
can form a charmed baryon. Therefore, we may expect a small excess of over 
mesons. However, there are other subtle QCD effects that might contribute to a charm 
meson production asymmetry |3l H] ; we note for b quarks the asymmetries are estimated 
to be at the 1% level [5], and we would expect them to be smaller for c quarks, although 
quantitative predictions are difficult. Another conjecture is that any asymmetries might 
be reduced as particles are produced at more central rapidities. 




Figure 1: Production of cc quark pairs in a pp collision via gluons. 

Measurements of CP violating asymmetries in charm and bottom decays are of prime 
importance. These can be determined at the LHC if production and detection asym- 
metries are known. The measurement of asymmetries in flavour specific modes usually 
involves detection of charged hadrons, and thus requires the relative detection efficiencies 
of TT"'" versus 7i~ or versus K~ to be determined. While certain asymmetry differ- 
ences can be determined by cancelling the detector response differences to positively and 
negatively charged hadrons [6], more CP violating modes can be measured if the relative 
detection efficiencies can be determined. 

In this Letter we measure the production asymmetry. 



1 



where (j{Dj) is the inclusive prompt production cross-section. We use Df — )■ (pn"^ decays, 
where (p ~^ K^K~ . Since Df — )■ 07r^ is Cabibbo favoured, no significant CP asymmetry 
is expected [71 E] . Assuming it to be vanishing, A-p is determined after correcting for the 
relative Df and D~ detection efficiencies. Since the final states are symmetric in kaon 
production, this requires only knowledge of the relative vr"*" and vr" detection efficiencies, 
e(vr+)/e(vr-). 

2 Data sample and detector 

The data sample is obtained from 1.0 fb~^ of integrated luminosity, collected with the 
LHCb detector [H] using pp collisions at a center-of-mass energy of 7 TeV. The detector 
is a single-arm forward spectrometer covering the pseudorapidity range 2 < 77 < 5, de- 
signed for the study of particles containing 6 or c quarks. The detector includes a high 
precision tracking system consisting of a silicon-strip vertex detector surrounding the pp 
interaction region, a large-area silicon-strip detector located upstream of a dipole magnet 
with a bending power of about 4Tm, and three stations of silicon-strip detectors and 
straw drift-tubes placed downstream. The combined tracking system has a momentum 
resolution Ap/p that varies from 0.4% at 5GeV to 0.6% at 100 GeVQ Charged hadrons 
are identified using two ring-imaging Cherenkov (RICH) detectors. Photon, electron and 
hadron candidates are identified by a calorimeter system consisting of scintillating-pad and 
pre-shower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons 
are identified by a muon system composed of alternating layers of iron and multiwire 
proportional chambers. The trigger consists of a hardware stage, based on information 
from the calorimeter and muon systems, followed by a software stage which applies a full 
event reconstruction. Approximately 40% of the data was taken with the magnetic field 
directed away from the Earth (up) and the rest down. We exploit the fact that certain 
detection asymmetries cancel if data from different magnet polarities are combined. 

Events are triggered by the presence of a charm hadron decay. The hardware trigger 
requires at least one hadronic transverse energy deposit of approximately 3 GeV. Subse- 
quent software triggers and selection criteria require a subset of tracks to not point to a 
primary pp collision vertex (PV), and form a common vertex. 

3 Measurement of relative pion detection efficiency 

In order to measure e(7r+)/e(7r~), we use the decay sequence D*^ — >■ vr+Z^^, D^ — )■ 
K'n^-K^n' , and its charge-conjugate decay, where 7r+ indicates the "slow" pion com- 
ing directly from the D*^ decay. Assuming that the D*^ comes from the PV, there are 
sufficient kinematic constraints to detect this decay even if one pion from the D^ decay is 
missed. We call these "partially" reconstructed decays. We can also "fully" reconstruct 
this decay. The ratio of fully to partially reconstructed decays provides a measurement of 

"'^We work in units with c—1. 
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the pion reconstruction efficiency. We examine D*~^ and D*^ candidate decays separately, 
and magnet up data separately from magnet down data. The latter is done to test for 
any possible left-right detector asymmetries. In both cases the missing pion's charge is 
required to be opposite of that of the detected kaon. 

Kaon and pion candidates from candidate decays are required to have transverse 
momentum, > 400 MeV, and a track quality fit with P^r number of degrees of 
freedom (ndf) < 3, keeping more than 99% of the good tracks. The distance of closest 
approach of track candidates to the PV is called the impact parameter (IP). A restrictive 
requirement is imposed on the IP x^? which measures whether the track is consistent with 
coming from the PV, to be greater than 4. In addition both particles must be identified in 
the RICH. For the vr^, the px requirement is lowered to 250 MeV, with both IP < 0.3 mm 
and IP < 4 being required. Further tight restrictions are placed on candidates. 
The candidate tracks from the decay must fit to a common vertex with x^/ndf < 6, the 

candidates must have a flight distance of at least 4 mm from the PV and have a flight 
distance > 120. We require 1.4 < m{K~n~^7r~) < 1.7 GeV, and that the invariant 
mass of the it~^tt~ candidates must be within ±200 MeV of the p(770) mass, to improve 
the signal to background ratio. 

We select partially reconstructed right-sign (RS) candidates by examining the 
mass difference Amprt = m{TT'^K^'K^7i^) — m{K^'K^'K^). Wrong-sign (WS) candidates 
are similarly selected but by requiring that the charge of the kaon be the same as that of 
the 7r+. Figure [2] shows distributions of Amprt for magnet up data. Note that the yield 
of WS events is reduced due to a prescale factor applied in the selection. 




Figure 2: Distributions of mass differences in partial reconstruction for (a) RS 
m{'K^ K~ T{^) — m{K^7T^'K^) and (b) WS m{7T'^K~^7r~^7T~) — m(i^"'"7r"'"7r~) candidates, 
for magnet up data. The (green) dotted line shows the signal, the (red) dashed line the 
background, and the (blue) solid line the total. The fit shapes are defined in Appendix 
A. 
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In order to determine the size of the signals above the background we perform simulta- 
neous binned maximum likelihood fits to the RS and WS distributions. The parametriza- 
tion of the signal probability density function (PDF) is given in Appendix A. The signal 
and background PDFs are identical for RS and WS D° and D° events, only the absolute 
normalizations are allowed to differ. We also include a "signal" term in the fit to WS 
events to account for the doubly-Cabibbo-suppressed (DCS) signals. The ratio of the 
DCS signal in WS events to the signal in RS events is fixed to that obtained in the mass 
difference fit in full reconstruction. 

Using momentum and energy conservation and knowledge of the direction of the 
flight direction, the inferred three-momentum of the missing pion, Pi^i, is reconstructed us- 
ing a kinematic fitting technique [10] . Our resolution on inferred pions may be determined 
from the the fully reconstructed D*^ sample, by removing one detected pion whose three- 
momentum is well known, and treating the track combination as if it was partially recon- 
structed. We then have both detected and inferred momentum, and thus a measurement 
of the missing pion momentum resolution distribution, AP/P = (Pdetected — -Pinf)/-Pdetected• 
For further study, we take only combinations with good inferred resolution by accepting 
those where i^nf divided by its calculated uncertainty is greater than two and also where 
the transverse component of Pi^{ divided by its uncertainty is greater than 2.5; this elim- 
inates about 37% of the sample. 

In our sample of partially reconstructed events, we subsequently look for fully recon- 
structed decays by searching for the missing track. Candidate tracks must have p > 2 GeV, 
Pt > 300 MeV and be identified as a pion in the RICH. They also must form a vertex with 
the other three tracks from the decay with a vertex fit x^/iidf< 6, and have a four-track 
invariant mass within 30 MeV of the mass peak. Certain areas of the detector near its 
edges preferentially find only one charge or the other depending on magnet polarity. We 
remove fully reconstructed candidates where the detected pion projects to these regions, 
discarding 3% of the candidates. 

The mass difference for fully reconstructed combinations, Amfuu = 
■m(TTfK^-K^-K^-K^) — ■m(K^TT^-K^-K^) IS showu in Fig. [3| for both RS and WS cases. Only 
D*^ data in the magnet up configuration are shown. The shape of the mass difference 
signal PDF is described in Appendix B. 

In order to extract the signal yields, we perform a binned maximum likelihood 
fit to the D*^ and D*~ events, both RS and WS, simultaneously. Table [l] lists the 
signal yields for partial reconstruction, A^prt, and full reconstruction, iVfuu. The effi- 
ciency ratios are derived from ratios of RS yields, e{Tr~^) = A^fuii(-D°7r+)/A''prt(-D°7r+) and 
e{n~) = Nf^^\l{D'^'^^~) / Npj.t{D^7ij) . (The absolute efficiency inferred from these yields in- 
cludes geometric acceptance effects.) The p and p^ spectra of the tt^ used for the efficiency 
measurement are shown in Fig|4j 

The ratios of pion detection efficiencies are 0.9914±0.0040 and 1.0045±0.0034 for mag- 
net up and magnet down, respectively, with statistical uncertainties only. To obtain the 
efficiency ratio as a function of momentum we need to use the inferred momentum of the 
missing pion. Because of finite resolution it needs to be corrected. This is accomplished 
through an unfolding matrix estimated using the fully reconstructed sample by comparing 
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Figure 3: Distributions of the mass difference Amfuu for (a) RS and (b) WS events using 
magnet up data. The (green) dotted hne shows the signal, the (red) dashed hne the 
background, and the (blue) solid line the total. The fit shapes are defined in Appendix 
B. 

Table 1: Event yields for partial and full reconstruction. 

Category Magnet up Magnet down 

A^prt(^°vr+) 460,005±890 671,638±1020 

Np,t{D^7v;) 481,823±873 694,268±1035 

Ni^nlD%+) 207,504±465 299,629± 570 

Nf^ir{D^Tt-) 219,230±478 308,344± 579 



the measured momentum of a found pion that is then ignored and its momentum inferred 
using the kinematic fit. The efficiency ratio is shown as a function of momentum in Fig. |5j 
Most systematic uncertainties cancel in the efficiency ratio, however, some small residual 
effects remain. To assess them we change the signal and background PDFs in full and 
partial reconstruction by eliminating, in turn, each of the small correction terms to the 
main functions. The full fit is then repeated. Each change in the efficiency ratios are 
between 0.01—0.02%. We also change the amount of DCS decays by the measured uncer- 
tainty in the branching fraction. This also gives a 0.020% change. The total systematic 
error is 0.045%. Furthermore, the entire procedure was checked using simulation. 

Although we correct relative pion efficiencies as a function of p, it is possible that 
there also is a pt dependence that would have an effect if the pt distributions of the D*"^ 
and Df were different. The efficiency ratios for different slices are shown in Fig. [6] For a 
fixed p interval there is no visible pi dependence. 
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Figure 4: Distribution of fully reconstructed signal candidates for magnet up data as a 
function of pion (a) p and (b) px- 
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Figure 5: Relative detection efficiency in bins of detected pion momentum: (red) circles 
represent data taken with magnet polarity up and (blue) squares show data taken with 
magnet polarity down. Only statistical errors are shown. 

The relative pion efficiencies are consistent with being independent of p and px- The 
tracking acceptance does depend, however, on the azimuthal production angle of the 
particles, ip. This is mostly because tracks can be swept into the beam pipe and not be 
detected by the downstream tracking system. Therefore, for purposes of the production 
asymmetry analysis we determine e(7r+)/e(7r^) as a function of in two momentum 
intervals: 2 — 20 GeV, and above 20 GeV. The r.m.s. resolution on the inferred ip is 
0.25 rad, much smaller than the 7r/4 bin size. The correction factors are shown in Fig. [7j 
The average correction for magnet up and magnet down is consistent with unity. Thus 
any residual biases in the Df yields due to n^/n~ asymmetries will also cancel in the 
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Figure 6: Relative efficiency averaged over magnet up and magnet down samples versus 
pion p and px- 
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Figure 7: Azimuthal angle distribution of e(7r+)/e(7r^) for magnet up data (red circles) 
and magnet down data (blue squares), and their average (black diamonds) for (a) pion 
momentum 2 < p < 20 GeV and (b) p > 20 GeV. 



4 production asymmetry 

The decay Df — )• K'^K'n^ is used with the invariant mass of the K^K^ required to be 
within ±20 MeV of the mass. Events are triggered at the hardware level by requiring 
that either the or the K~ deposits more than 3 GeV of transverse energy in the hadron 
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Figure 8: Invariant mass distributions for (a) K^K^tt^ and (b) K^K^tx^ candidates, 
when m{K^K^) is within ±20 MeV of the mass, for the dataset taken with magnet 
polarity down. The shaded areas represent signal, the dashed line the background and 
the solid curve the total. 



calorimeter. Subsequent software triggers are required to select both decay products. 

To select a relatively pure sample of Df — )■ K~^K~7i'^ candidates each track is required 
to have xV^df <4, > 300 MeV, IP > 4, and be identified in the RICH. All three 
candidate tracks from the Df have pt > 2 GeV, and must form a common vertex that 
is detached from the PV. The requiring all three tracks to come from a common 
origin must be < 8.33, this decay point must be at least 100 /im from the PV, and 
the significance of the detachment must be at least 10 standard deviations. The Df 
candidates' momentum vector must also point to the PV, which reduces contamination 
from 6-hadron decays to the few percent level. We remove signal candidates with pions 
which pass through the detector areas with large inherent asymmetries, as we did to 
measure the relative pion efficiencies. 

Figure [s] shows the invariant mass distributions for (a) K'^K^n^ and (b) K^K~tc~ 
candidates for data taken with magnet polarity down. We perform a binned maximum 
likelihood fit to extract the signal yields. The fitting functions for both and Df 
signals are triple Gaussians where all parameters are allowed to vary, except two of the 
Gaussians are required to have the same mean. The background function is a second 
order polynomial. The numbers of Df events obtained from the fits are listed in Table |2} 

The rapidity of the Df is defined as 



1, E + p, 
?/ = ;r In 



E 



Pz 



(2) 



where E and Pz are the energy and z component of the Df momentum. We measure 
the production asymmetry Ap as a function of both Df y and Pt- In each y or p^ bin 
we extract the efficiency corrected ratio of yields by applying corrections as a function of 
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y Pt (GeV) 

Figure 9: (a) Df rapidity distribution (b) Df pt distribution for background subtracted 
magnet up data. The statistical uncertainty on the number of events in each bin is smaller 
than the line thickness. 



azimuthal angle in the two pion momentum intervals defined previously. Magnet up and 
down data are treated separately. The y and distributions are shown in Fig. [9j Here 
sidebands in KKtt mass have been used to subtract the background, where the sidebands 
are defined as between 30 — 70 MeV above and below the peak mass value of 1969 MeV. 
As this interval is twice as wide as the signal peak, we weight these events by a factor of 
1/2. 



Figure 10 shows Ap as a function of either y or px- The error bars reflect only the 
statistical uncertainties, which includes both the statistical errors on e{n~^) / e{7i~) and the 
Df yields; the error bars are partially correlated, the uncertainties from the Df yields 
are about half the size of those shown. The values in and y intervals are listed in 
Table [3j An average asymmetry in this y and p^ region can be derived by weighting 
the asymmetry in each bin by the production yields. Thus we take the asymmetry in 
each y and pt interval, weight by the measured event yields divided by the reconstruction 
efficiencies. The resulting integrated production asymmetry Ap is (—0.20 ± 0.34)%, and 
(—0.45 ± 0.28)%, for magnet up and magnet down samples, respectively. The errors are 
statistical only. Averaging the two results, giving equal weight to each to cancel any 
residual systematic biases, gives 

Ap = (-0.33 ± 0.13 ± 0.18 ± 0.10)%, 



Table 2: Fitted numbers of Df events for both magnet up and down data. 





Magnet up 


Magnet down 


Df 


152,696±448 


230,860±514 


D- 


154,209±438 


233,266±549 
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Figure 10: Observed production asymmetry Ap as a function of (a) y, and (b) p^. The 
errors shown are statistical only. 



where the first uncertainty is statistical from the Df yields, the second statistical due to 
the error on the efficiency ratio and the third systematic. The systematic uncertainty on 
Ap has several contributions. Uncertainties due the background shape in the Df mass 
fit are evaluated using a higher order polynomial function, that gives a 0.06% change. 
Statistical uncertainty on MC efficiency adds 0.06%. Constraining the signal shapes of 
the Df and D~ to be the same makes a 0.04% difference. Possible changes in detector 
acceptance during magnet up and magnet down data taking periods are estimated to 
contribute 0.03%. The systematic uncertainty from the pion efficiency ratio contributes 
0.02%. Differences in the momentum distributions of K~ and that arise from inter- 
ference with an S-wave component under the (p peak can introduce a false asymmetry 
[TT] . For our relatively high momentum Df mesons this is a 0.02% effect. Contamination 
from b decays causes a negligible effect. Adding all sources in quadrature, the overall 
systematic uncertainty on Ap is estimated to be 0.10%. 



Table 3: Ap (%) shown as a function of both y and px- 



Pt (GeV) 




y 






2.0-3.0 


3.0-3.5 


3.5-4.5 


2.0- 6.5 


0.2 ±0.5 


-0.7 ±0.5 


-0.9 ±0.4 


6.5- 8.5 


-0.3 ±0.4 


0.1 ±0.5 


-1.2 ±0.5 


8.5-25.0 


0.2 ±0.3 


-0.3 ±0.5 


-1.0 ±0.8 
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5 Conclusions 



We have developed a method using partially and fully reconstructed D*"^ decays to mea- 
sure the relative detection efficiencies of positively and negatively charged pions as a 
function of momentum. Applying this method to Df mesons produced directly in pp 
collisions, i.e. not including those from decays of h hadrons, we measure the overall 
production asymmetry in the rapidity region 2.0 to 4.5, and p^ > 2 GeV as 

= ""i nl? 7 ""i \ = (-0-33 ± 0.22 ± 0.10)%. (3) 

The asymmetry is consistent with being independent of px, and also consistent with being 
independent of y, although there is a trend towards smaller A-p values at more central 
rapidity. These measurements are consistent with theoretical expectations [31 H], provide 
significant constraints on models of Df production, and can be used as input for CP 
violation measurements. 



Acknowledgments 

We express our gratitude to our colleagues in the CERN accelerator departments for 
the excellent performance of the LHC We thank the technical and administrative staff at 
CERN and at the LHCb institutes, and acknowledge support from the National Agencies: 
CAPES, CNPq, FAFERJ and FINEP (Brazil); CERN; NSFC (China); CNRS/IN2P3 
(France); BMBF, DEC, HGF and MPG (Germany); SFI (Ireland); INFN (Italy); FOM 
and NWO (The Netherlands); SCSR (Poland); ANCS (Romania); MinES of Russia and 
Rosatom (Russia); MICINN, XuntaGal and GENCAT (Spain); SNSF and SER (Switzer- 
land); NAS Ukraine (Ukraine); STFC (United Kingdom); NSF (USA). We also acknowl- 
edge the support received from the ERG under FP7 and the Region Auvergne. 



Appendix A: Fitting functions for partial reconstruc- 
tion 

The signal probability density function (PDF) is given by: 

/sig(Amprt) = /eff(Amprt) • BG{Am^,t\ cru o^r, /i), where (4) 
^^(Amprt; /i, ai) if Amp^t < /i, 



SG(Amprt) 



r^G(Amprt; /i, a^) if Amprt > yW- 



G(Amprt;/i,o") is a Gaussian function with mean /i and width a, and BG{/S.m^,-t) is a 
bifurcated Gaussian function. The efficiency function /eff(Amprt) is defined as: 



|a(Amprt-Amo)|^ 



if Amprt ~ Amo > 0, 



/eff(Amprt) = <^ l+l-{Amp..-A„^o)r ^"^P^^* " (5) 

if Amprt - Amo < 0, 
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where a, N, and Amo are fit parameters. The resolution function (the bifurcated Gaussian 
function) is multiphed by the efficiency function /egf(Amprt) in order to account for the 
"turn-ofF' behaviour of the quantity Amprt near the threshold (pion mass) . There are in 
total six shape parameters in this signal PDF which arc left to vary in the fit. 

The background PDF is taken as a threshold function with the inclusion of extra 
components to obtain a good description of the WS combinations. It is defined similarly: 

/bkg(Amprt) = r(Amprt)-(c2Am2rt + ciAmprt + l)-/i-SG'(Amprt) + /2-G'(Amprt), (6) 

/*(Amprt) = [1 - exp(-(Amprt - Am^)/cp)] ■ a^--*/^< + bpiArripJAml - 1). (7) 

The parameters used in the background functions BG and G are different than the ones 
used in the signal functions. There are in total 11 shape parameters in the background 
PDF that are determined by the fit. We also fit using /* (Amprt) as the background PDF 
alone to estimate the systematic uncertainty on the efficiency ratio. 



Appendix B: Fitting functions for full reconstruction 

The signal PDF is defined as: 

/sig(Amfuii) = /iG'(Amfuii; //i, ai) + /2G'(Amfuii; //2, (72) 

+ (1 - /l - /2)/student(Amfuli; Awo, Ul, Uh: CT^^ve: Sa), (8) 

where G(Amfuii) is a Gaussian function defined in Appendix A, and /student (A'^fuii) is 
obtained from the Student's t-distribution 

_ r(./2 + 1/2) / tn(-'^/-V2) 

where F is the Gamma function. We define t = (Amfuu — Amo)/o' with Amo and a the 
mean and width. In order to obtain the asymmetric t-function, the width parameter a 
and number of degrees of freedom i> are allowed to be different for the high and low sides 
of Amfuu. Widths for high and low sides of Amfuu are then defined as: ah = o"ave + ^o-, 
and ai — a^ve — Sa, and u parameters for high and low sides are denoted as i/h and ui, 
respectively. The bifurcated Student's t-function can then be defined as: 




/student (Amfull) = { ^ .a™,..„-a™o)2\ {->^il2-l/2) (10) 



if Amfull — Amo > 0, 
if Amfull — Amo < 0. 



Auxiliary terms arc defined as: 



r(z/j2 + l/2) r(z/,/2 + l/2) 2pi 2ph 

Pf" ^ r? — aTV^ — I ' = p/ /ON ^1 P = — — , ri = — — . (11) 
T{uh/2)^/i7i,\ah\ T{ui/2)^i\ai\ Ph + Pi Ph + Pi 
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In total there are 11 shape parameters in the signal PDF, all of them are allowed to vary 
in the fit. The background PDF is extracted from WS events, and is defined as: 

/bkg(AmMi) = (1 - h) ■ /*(AmMi) + /s ■ fiG'(AmMi), (12) 

where /*(Amfuii) is defined in Eq. [7| We add a correction function, a bifurcated Gaussian, 
in order to have a better fit; the shape and the fraction of the bifurcated Gaussian is 
determined empirically from WS events. (We also use a background shape without this 
correction term to estimate the systematic uncertainty on the efficiency ratio.) We include 
a "signal" term in the fit to WS events to account for doubly-Cabibbo suppressed decays. 
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